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Abstract 

Variations in elevation limit tlie growtli and distribution of alpine plants because multiple environmental stresses impact 
plant growth, including sharp temperature shifts, strong ultraviolet radiation exposure, low oxygen content, etc. Alpine 
plants have developed special strategies to help survive the harsh environments of high mountains, but the internal 
mechanisms remain undefined. Kobresia pygmaea, the dominant species of alpine meadows, is widely distributed in the 
Southeastern Tibet Plateau, Tibet Autonomous Region, China. In this study, we mainly used comparative proteomics 
analyses to investigate the dynamic protein patterns for K. pygmaea located at four different elevations (4600, 4800, 4950 
and 5100 m). A total of 58 differentially expressed proteins were successfully detected and functionally characterized. The 
proteins were divided into various functional categories, including material and energy metabolism, protein synthesis and 
degradation, redox process, defense response, photosynthesis, and protein kinase. Our study confirmed that increasing 
levels of antioxidant and heat shock proteins and the accumulation of primary metabolites, such as proline and abscisic 
acid, conferred K. pygmaea with tolerance to the alpine environment. In addition, the various methods K. pygmaea used to 
regulate material and energy metabolism played important roles in the development of tolerance to environmental stress. 
Our results also showed that the way in which K. pygmaea mediated stomatal characteristics and photosynthetic pigments 
constitutes an enhanced adaptation to alpine environmental stress. According to these findings, we concluded that K. 
pygmaea adapted to the high-elevation environment on the Tibetan Plateau by aggressively accumulating abiotic stress- 
related metabolites and proteins and by the various life events mediated by proteins. Based on the species'lexible 
physiological and biochemical processes, we surmised that environment change has only a slight impact on K. pygmaea 
except for possible impacts to populations on vulnerable edges of the species' range. 



•0-PLOS I o-^E 



Citation: Li X, Yang Y, Ma L, Sun X, Yang S, et al. (2014) Comparative Proteomics Analyses of Kobresia pygmaea Adaptation to Environment along an Elevational 
Gradient on the Central Tibetan Plateau. PLoS ONE 9(6): e98410. doi:10.1371/journal.pone.0098410 

Editor: Eryuan Liang, Chinese Academy of Sciences, China 

Received January 13, 2014; Accepted May 1, 2014; Published June 2, 2014 

Copyright: © 2014 LI et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the Major State Basic Research Development Program (201 0CB95 1700), the National Natural Sciences Foundation of China 
(31 170256) and the Young Academic and Technical Leader Raising Foundation of Yunnan Province (2012HB041). The funders had no role in study design, data 
collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: huxiangyang@mail.klb.ac.cn (XH); yangyp@mail.klb.ac.cn (YPY) 



Introduction 

High elevation areas have always attracted the attention of 
ecologists, especially in light of global climate change. The unique 
alpine environments affect biological survival and evolution [1] 
and high-elevation regions are considered to be more sensitive to 
climate change than other areas [2,3]. Most alpine plants 
experience the harsh climates in the high mountain elevations, 
such as extremely low temperatures, low oxygen concentrations, 
strong ultraviolet (UV) radiation, aridity and frequendy violent 
wind conditions. Therefore, alpine plants have evolved with many 
defense mechanisms designed to help them to survive at high 
elevations. For example, some alpine plants exhibit dwarfism with 
a small leaf area index, which avoids stress caused by strong wind 
[4] . A hairy floss covers the leaves and stems of most alpine plants 
to protect them from low temperatures [5]. Plants can also avoid 
water loss by deepening root structures and increasing stomatal 



control [5]. Previous studies demonstrate that alpine plants use 
several morphological strategies to adapt to the high-elevation 
environment, but the underlying proteome and physiological 
mechanisms of the adaptation remain undefined [6] . In addition, 
physiological and biochemical characteristics are more sensitive to 
environmental change when compared with more stable mor- 
phology characteristics. Therefore, studies in these directions are 
essential because the impact of climate change on alpine plants 
may be predictable. 

Kobresia pygmaea is a dominant species of alpine meadows in the 
southeastern humid Tibetan highlands and the high alpine 
pastures of the southern and eastern declivity [7]. Widespread 
along the alpine elevation gradient, it is found from the lowest 
slopes (3000 m) to the highest outposts (5960 m) [7], indicating its 
successful acclimation to the alpine environment. K. pygmaea has 
many excellent features making it viable at many microsites, 
including resistance to low temperature, drought, trampling, and 
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soil erosion, which are all of great importance in maintaining the 
stability of the regional environment [8] . Therefore, it is greatly 
important to understand interactions between K. pygmaea and 
environments. Elevation gradient causes gradient change of 
environment factors and is considered to be the epitome of 
horizontal zone. This supplies a good chance for studying alpine 
plants adaptation to environment conditions. We hypothesized 
that environmental conditions became more severe with elevation 
rising; and what physiological and proteomics characteristics K. 
pygmaea applied to survive? In the present study, we selected four 
elevational gradients to understand the integral adaptation 
mechanisms of K. pygmaea. The results showed that K. pygmaea 
used multiple strategies to acclimate itself to alpine environmental 
stress; precisely because of this, we hypothesized environment 
change might have litde impact on K. pygmaea across a wide area. 

Results 

Vertical investigation and analysis of environmental 
factors 

In our sampUng sites (30°30'-30°32'N, 91°03'E), the alpine 
plant K. pygmaea is mainly distributed from 4600 to 5200 m in 
elevation; however, as the elevation reaches .5200 m the popula- 
tions fall to a very low level. Various factors, including soil 
temperature, humidity, average sunlight time, average rainfall and 
average UV radiation exposure, vary along the elevation gradient. 
Ma et al. [5] found soil temperatures decreased significantly (P< 
0.05) and rainfall increased with increasing elevation at this 
location during July. Principal component analysis indicated that 
soil temperature, precipitation (or soil water content) and average 
UV radiation were three main factors that limited the growth and 
distribution of K. pygmea. Therefore, we mainly considered the 
effects of temperature, precipitation and solar radiation on K. 
pygmaea adaptation along an elevation gradient in our experiments. 

Variation of stomatal density and stomatal aperture 
length with increasing elevations 

The stomatal density of K. pygmaea increases gradually with 
elevation from 4600 to 4950 m reaching a maximum at 4950 m 
and then decreasing at 5 1 00 m to a level similar to that at 4800 m 
(Figure 1 A left). In contrast, the stomatal aperture length obviously 
declined along the elevation gradient, despite littie difference 
observed between 4800 and 4950 m (Figure lA right & B). 

Differences of photosynthetic pigments content along 
an elevational gradient 

The total chlorophyll content of K. pygmaea rose significantiy 
from 4600 to 4950 m; it increased only slightly from 4950 to 
5100 m (Figure 2A left). Carotenoid content rose continuously 
with increasing elevation from 4600 to 5100 m (Figure 2A right). 

Changes of proline and abscisic acid content with 
increasing elevation 

The proline and abscisic acid (ABA) content both rose obviously 
with increasing elevation from 4600 to 5100 m (Figure 2B). 

Dynamics of the antioxidant enzyme system along the 
elevation gradient 

To investigate reactive oxygen species (ROS) metabolism of K. 
pygmaea along the elevational gradient, we detected catalase (CAT; 
EC 1.11.1.6), ascorbate peroxidase (APX; EC 1.11.1.11), gluta- 
thione reductase (OR; EC 1.8.1.7), and superoxide dismutase 
(SOD; EC 1.15.1.1) activities. The results revealed that they 



exhibited analogous variation. They all increased from 4600 to 
5100 m, although the amplitude of change varied at different 
elevations (Figure 3). 

Protein profiling of K. pygmaea adaptation to the alpine 
environment with increasing elevation 

To investigate the response of proteins in K. pygmaea to 
elevational change, we performed two-dimensional electrophoresis 
(2-DE) to identily the entire protein accumulation profile in K. 
pygmaea from 4600 to 5100 m. We performed three biological 
replicates and gels were visualized by coomassie brilliant blue 
(CBB) staining (Figure 4, Figure SI and S2 in File SI). After 
staining, proteins were analyzed by PDQuest software (Bio-Rad, 
Hercules, CA, USA). Under stringent conditions, all differentially 
displayed proteins were unambiguously identified by MALDI- 
TOF-MS/MS analysis and compared to the NCBI non-redun- 
dant database. 58 protein spots were successfully quantified with 
significant differential expression changes among the leaf samples 
from four different elevation locations (4600, 4800, 4950 and 
5100 m). All of the proteins demonstrated increased expression (> 
1.50) or decreased expression (<0.67) at high elevations (4800, 
4950 and 5100 m) compared to samples from the lowest elevation 
(4600 m), which was used as the control (Table 1 and Figure 4). 
According to NCBI gene annotations, the identified proteins could 
be classified into six functional groups: material and energy 
metabolism, protc'iii synthesis and degradation, redox process, 
photosynthesis, defense response, and protein kinase (Figure 5A 
and Table 1). A hierarchical cluster analysis was conducted to 
categorize the proteins that showed differential expression profiles 
in response to the elevational gradient (Figure 5B). More 
obviously, antioxidant enzymes, defense response proteins, pho- 
tosynthesis-related proteins and some material and energy- 
associated proteins were observed to cluster together, which were 
up-regulated with increasing elevation (Figure 5B and Table 1). 
Venn diagram analysis was used to reflect different change 
patterns of proteins at each higher elevation compared with 
4600 m, respectively (Figure 5C), which further proved that many 
more proteins in K. pygmaea were up-regulated rather than down- 
regulated with an increase in elevation. 

Changes of malondialdehyde content with increasing 
elevation 

The content of malondialdehyde (MDA) changed slightiy and 
insignificantiy from 4600 to 5100 m, (Figure 6A). 

Differences in the expression of specific proteins with 
elevational change 

The results of western blot analysis showed that 9-cis- 
epoxycarotenoid dioxygenase (NCED), an enzyme regulating 
ABA synthesis, had continuously increasing expression from 
4600 to 5100 m (Figure 6B). Similar variations were found in 
the expression of dehydrin and heat shock protein 18.2 (HSP18.2) 
(Figure 6B). The expression of mitogen-activated protein kinase 6 
(MAPK6) appeared to be different from these; it increased to a 
peak from 4600 to 4950 m but fell at the elevation of 5100 m 
(Figure 6B). 

Discussion 

Coordination between stomatal size and density at high 
elevations 

Stomata allow plants to exchange gases with the surrounding 
atmosphere. Stomatal density and aperture directiy affect the 
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Figure 1 . Stomatal characteristics of K. pygmaea from different elevations. A: Changes in stomatal density (left) and aperture length (right) 
In K. pygmaea along an elevational gradient. B: Stomatal shape and size changes of K. pygmaea at different elevations at the same magnification. 
Error bars indicate SB. Means denoted by different letters are significantly different (P<0.05) (A). 
doi:l 0.1 371 /journal.pone.009841 0.gOOl 



photosynthetic and transpiration rates of plants [9]. Most 
environment factors can affect the growth and conductance of 
stomata [10—12]; but the relationships between environment 
factors and stomatal traits are still poorly known bec:ause of their 
multiple correlations in different species [11—14]. According to 
past research [15-18], no exact relationship exists between 
stomatal characteristics and elevational gradients in various 
species. 

In the present study, temperature decreased and moisture 
increased from 4600 to 5100 m, while stomatal aperture shrank 
and stomatal density initially rose (4600-4950 m) but then 
declined (5100 m) (Figure lA). One reasonable explanation for 
reduced stomatal aperture is that smaller stomatal sizes are 
beneficial for K. pygmaea, allowing the plant to avoid damage from 
lower temperature and stronger UV radiation with increasing 
elevation. However, at the lower altitude, because K. pygmaea needs 
have some degree of drought tolerance due to the lower rates of 
precipitation, larger stomatal aperture does not help to reduce 
water loss; thus, K. pygmaea remedies this problem by decreasing 
stomatal density. In addition, environmental conditions may be 
more unstable at higher elevations, implying a higher frequency of 
stomatal opening and closing. Because K. pygmaea leaves can curl 
into a needle-shape in the wild, we presume smaller stomatal 
apparatuses are more flexible in situations where plants have 
limited leaf areas. That is why the stomatal aperture length 
significandy decreased at 5100 m (Figure lA left); nevertheless, to 
maintain adequate gas exchange space, stomatal density increases, 
despite a decrease at 5100 m. Speculation may provide some 
indirect evidence from the accumulation of ABA (Figure 2B left) 
and its synthase NCED (Figure 6B), since ABA plays an important 



role in stomatal opening and closing [19,20]. MAPK signal 
pathway has also been reported to be involved in the regulation of 
stomatal development and movement [21,22]. Interestingly, the 
expression of MAPK6 was observed to be consistent with variation 
of stomatal density (Figure 6B), indicating their close relationship. 

Most importandy, stomatal development must be affected both 
positively and negatively by different environmental factors. One 
result of adaptation to the environment is the coordination 
between stomatal size and density. 

Photosynthetic pigments and solar radiation 

A plant's photosystem is the basis of photosynthesis, and 
photosynthetic pigments are mainly involved in light absorption 
and transmission. As elevation increases, solar radiation is 
enhanced. While this supplies plants with more light, excessive 
Kght may damage the photosystem of plants [23]. In our results, 
two main photosynthetic pigments, chlorophyll and carotenoid, 
both showed an increasing trend from 4600 to 5100 m (Figure 2A); 
interestingly, several proteins related to light capture and 
absorption (spots 36, 37, 39, 54 and 58) displayed similar 
expression variation (Figure 5B). Chlorophyll content increased 
significandy from 4600 to 4950 m but only slighdy at higher 
elevations (Figure 2A left). This implies chlorophyll could absorb 
enough hght for photosynthesis; however, excessive hght must be 
prevented to avoid damaging R. pygmaea. This may he llic reason 
that chlorophyll content remained steady between 4950 and 
5100 m. Carotenoids, as a class of accessory pigments, showed 
sustained increase from 4600 to 5100 m (Figure 2A right). Apart 
Irom the normal function for photosynthesis, carotenoids are 
believed to help plants absorb excessive light to prevent 
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Figure 2. Changes of metabolite content in K. pygmaea with elevation increasing. A: Changes of photosynthetic pigment content in K. 
pygmaea at different elevations. B: Changes of abscisic acid and proline content in K. pygmaea along an elevational gradient. Error bars indicate SE. 
Means denoted by different letters were significantly different (P<0.05). 
doi:l 0.1 371/journal.pone.009841 0.g002 
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Figure 3. Changes of antioxidant enzyme activities in K. pygmaea a\ong an elevational gradient. Error bars indicate SE. Means denoted by 
different letters were significantly different (P<0.05). 
doi:1 0.1 371/journal.pone.009841 0.g003 
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Figure 4. Dynamic protein changes in leaves of K. pygmaea a\ong an elevational gradient. All 58 Identified protein spots were shown on 
the 2D gel at 4600 m. Other important protein spots specially mentioned appeared on the 2D gels of 4800, 4950 and 5100 m. 
doi:1 0.1 371 /journal.pone.009841 0.g004 



chlorophyll photo-oxidation by directly assimilating UV radiation 
in a way that limits damage to plants at high elevations [24,25]. 

The important role of the antioxidant enzyme system in 
K. pygmaea 

Reactive oxygen species (ROS) wiU be generated when plants 
undergo aerobic metabolism, e.g., photosynthesis and respiration 
[26]. If superabundant ROS are not removed promptly, plants 
may experience oxidative stress that may eventually lead cell to 
death [26]. Thus, plants have evolved scavenging machineries 
including antioxidant enzymes and antioxidants to keep ROS at 
physiologically acceptable levels [27]. 

Here, we determined the activities of four antioxidant enzymes 
(CAT, APX, SOD and GR) to investigate the response of K. 



pygmaea to the complex environment along an increasing 
elevational gradient. The results agreed with previous reports that 
some antioxidant enzyme activities were induced by low temper- 
ature, UV radiation or some other types of stresses [28-31]. We 
observed that all four enzyme activities exhibited similar increasing 
trends (Figure 3), which were consistent with the results related to 
another alpine herb, Potentilla saundersiana [5]. Interestingly, the 
results were identical for the expression of proteins in the 
antioxidant system (spots 42, 43, 46, 50 and 57) according to 
proteomics results (Figure 5B). On the contrary, MDA, which 
reflected grades of cellular oxidation [32], had litde change among 
the elevational gradient (Figure 6A). This fully certificated that no 
greater oxidative toxicity to K. pygmaea was caused by increasing 
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Figure 5. Comparative proteomics analyses results of four samples from different elevations. A: Functional classification of the identified 
proteins based on NCBI annotation. B: Hierarchical clustering of the identified protein expression profiles at different elevations. C: Venn diagram 
analysis of differentially expressed proteins at each higher elevations compared with 4600 m. Different colors correspond to the proteins' log- 
transformed fold-change ratios depicted in the bar at the bottom of the figure (B). 
doi:1 0.1 371 /journal.pone.009841 0.gOOS 



elevation, implying stimulated antioxidant enzymes played an 
important role in maintaining ROS homeostasis. 

Effects of material and energy-associated proteins 

Plants growing at higher elevations require relatively more 
energy to respond to serious environmental stress [33]. We 
observed that the expression of most proteins related to material 
and energy metabolism changed dramatically at high elevations. 
Carbohydrates provide the main nutrients used to supply organs 
with energy. In our experiments, the proteins involved in reducing 
or non-reducing sugar synthesis were difiFerentiaUy up-regulated 
with increasing elevation, including glyceraldehyde-3-phosphate 
dehydrogenase B subunit (spot 16), glyceraldehyde 3-phosphate 
dehydrogenase A subunit (spot 1 7), fructose-bisphosphate aldolase 
(spots 18 and 21), and sucrose phosphate phosphatase (spot 38). 
These findings agreed with the expression of ATP synthases (spots 
5 and 53) (Figure 5B), which regulate ATP synthesis to provide 
energy for plants. Meanwhile, ribulose bisphosphate carboxylase/ 
oxygenase activase, chloroplast precursor (spot 7), and triose 
phosphate isomerase cytosolic isoform (spot 44) showed differential 
down-regulation with increasing elevation (Figure 5B). In addition, 
other material and energy-associated proteins had irregular 
variation along the elevational gradient (Figure 5B), suggesting 
the regulation mechanisms used for carbon metabolism in K. 



pygmaea in response to alpine environmental stress are quite 
complicated. 

Glutamine synthetase (GS) is the key enzyme for plant ammonia 
assimilation, playing a significant role in normal life activities of 
plants. Higashi et al. [34] reported that lower GS activity would 
affect a variety of intracellular enzymes related to nitrogen 
metabolism and part of the glucose metabolism in radish. In the 
present study, glutamine synthetase la (spot 9) and glutamine 
synthetase (spots 12 and 13) were up-regulated with increasing 
elevation (Figure 5B), suggesting the significant regulation of 
glutamine for this alpine plant in response to environmental stress. 

The roles of abiotic stress-related metabolites and 
proteins 

Both proline and ABA can protect cells against chUl and other 
stressors at various stages of acclimation [35]. Proline helps plants 
avoid oxidative damage and indicates the presence of a stress 
response at the cellular level in many plants [36]. It is also believed 
to mediate osmotic adjustment, stabilize macromolecules, and 
store carbon and nitrogen for use during stress regimes [37]. ABA 
is an important plant hormone that modulates responses to abiotic 
stresses, including cold, heat and drought [38]. In our study, the 
content of both proline and ABA obviously increased when the 
elevation rose (Figure 2B). The environmental induction reflected 
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their important roles iii responding to low temperature, intense 
light and UV radiation, etc. with increasing elevation. 

ABA synthesis is partly regulated by NCED; and ABA pardy 
regulates the synthesis of dehydrins, which are also induced by 
cold, salt, and drought stress in plants [39] . Based on our western 
bolt results, the expression of NCED and dehydrin gradually 
increased at higher elevations, perfectly matching the change in 
ABA content (Figure 2B left & 6B). This not only confirmed their 
close relationships, but also verified the importance of NCED and 
dehydrin for the adaptation of K. pygmaea to the alpine 
environment. 

Heat shock proteins (HSPs) are an important group of 
protective proteins in eukaryotic organisms, which accumulate 
when plants are exposed to stressful conditions [40]. HSPs 
produced in plants can protect an organism's proteins from 
damage or repair damaged proteins, indicating that the induced 
formation of heat shock protein will help a plant acquire stress 
resistance [41]. In this study, both the results of proteomics (spots 2 
and 5 1) and western blot confirmed that rising elevations gradually 
increased the expression of HSPs in K. pygmaea (Figure 6B). This 
strongly suggested that a high level of HSP expression might 
confer K. pygmaea with greater tolerance to the complex 
environment of this alpine region. 

In addition, other defense response-related proteins, such as 
filamentation temperature-sensitive H 2B (spot 3), 14-3-3-like 
protein (spot 33) and NAC domain containing protein (spot 48), 
generally showed greater expression along the elevational gradient 
(Figure 5B). The findings supported the idea that the massive 
accumulation of abiotic stress-related metabolites and proteins 
improved the ability of K. pygmaea to adapt to the alpine 
environment. 

Effects of the protein synthesis and modification 

Proteins are the basis of life and metabolism, so the process of 
plant defense against harsh environments cannot be separated 
from the participation of various proteins. In our study, many 
proteins related to protein synthesis were detected with differential 
expression patterns along an increasing elevational gradient 
(Table 1 and Figure 5A). Protein post-translation modifications, 
like protein phosphorylation and dephosphorylation, can modu- 
late plant response to environmental stress [42] . Phosphoglycerate 
kinase, chloroplastic-like (spot 8), phosphoribulokinase (spot 1 1), 
nucleoside diphosphate kinase 1-like protein (spot 52), and 
nucleoside diphosphate kinase (spot 56) were differentially 
regulated with increasing elevation (Figure 5B), which demon- 
strated the complex regulation mechanisms K. pygmaea uses to 
respond to environmental stress. These findings suggested that 
protein synthesis and modification were closely related to the 
interaction between K. pygmaea and environment along increasing 
elevation. 

The predictable effects of environment change on K. 
pygmaea survival and development 

Environment is complex and changeable on the Tibetan 
Plateau, especially when climate change has become an important 
global issue that needs to be addressed; high-elevation areas are 
believed to be more sensitive to climate change than other areas 
[2,3] . Increasing temperature and changing precipitation patterns 
are two features of global climate change; they may interactively 
affect plant growth and development along an elevational 
gradient. 

In this study, temperature and precipitation show different 
combinations from the lowest to the highest primary range of K. 
pygmaea. Within its range of ecological amplitude, K. pygmaea 
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Figure 6. Changes of malondialdehyde content and special proteins expression in K. pygmaea with elevation increasing. A: Change of 
malondlaldehyde content in K. pygmaea along an elevational gradient. B: Protein accumulation of 9-cis-epoxycarotenoid dioxygenase, dehydrin, heat 
shock protein 18.2, and mitogen-activated protein kinase 6 in K. pygmaea from different elevations. Actin is included as a protein loading control. 
Error bars indicate SE. Means denoted by different letters were significantly different (P<0.05) (A). 
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successfully regulates its own physiological and biochemical 
processes to adapt to the environment, indicating that the species 
has a certain level of plasticity of adaptation to environmental 
changes. In another study, investigators analyzed how K. pygmaea 
successfully responded to the diurnal environment at an elevation 
of 4800 m at the protein level [43]. Taken together, whether on a 
spatial scale or time scale, K. pygmaea plants are able to manipulate 
the dynamic equilibrium of their internal environment to adapt to 
environmental change. Therefore, a certain extent of environ- 
mental change including climate change wiU not directly obstruct 
the survival and development of K. pygmaea, since climate change 
proceeds slowly; however, this idea does not address all issues 
related to this species. For example, K. pygmaea is limited by various 
environment factors on the fringe of its range, so perceptible 
environment change may disturb its ability to survive because the 
species is sensitive and vulnerable on the edge of its range. 
However, conclusive evidence of the species reaction to changing 
environment must be gained from long-term field observations or 
precise simulations of climate change. 

Conclusions 

Although K. pygmaea is subjected to different environmental 
conditions at various elevations across its range on the south-facing 
slope of the Nyainqentanglha Mountains, it is well acclimated to 
the alpine environment. We found several split-new strategies for 
K. pygmaea to adapt to environmental pressure resulting from the 
elevational gradient by a comparative proteomics and physiolog- 
ical approach. The strong accumulation of metabolites and 
proteins involved in abiotic stress and the many life events 
mediated by proteins were material foundations that help R. 
pygmaea acclimate to the alpine environment. Based on the 
flexibility of physiological and biochemical processes, we deemed 
that environment change would have a weak impact on K. pygmaea 
except for those vulnerable populations along the edge of the 
distribution. 

Complex and extensive interactions of various environment 
factors make it difficult to define the effect on K. pygmaea of single 
factor in the wild. This suggests that carefully controlled laboratory 
experiments would be necessary to further understand how K. 
pygmaea adapts to the environment. Nevertheless, natural environ- 



ments are not simple superpositions or offsets of the various factors 
that can be simulated in the lab; therefore, laboratory experiments 
cannot represent the actual status of this plant in the wild, even 
when the impact of each factor on plants is carefully analyzed and 
controlled. This fully demonstrates that field studies are also 
crucial to understanding the interaction between K. pygmaea and 
environment along an elevational gradient. 

Materials and Methods 

Ethics Statement 

This field study was permitted by the meadow managers. 
Institute of Geographic Sciences and Natural Resources Research, 
Chinese Academy of Sciences (CAS), Institute of Tibetan Plateau 
Research, CAS, and Grassland Station of Damxung County. The 
work was welcomed by Grassland Station of Damxung County 
and local herders. There wiU be no conffict of ethics and interest. 

Experimental samples collection 

This study was conducted in July 2012 on the south-facing slope 
of the Nyainqentanglha Mountains (30°30'-30°32'N, 9r03'E) 
near Damxung County, in central Tibet Autonomous Region, 
China (Figure 7 A and B). We selected one general location at an 
elevation ranging from 4300 to 5600 m where four sample sites 
were selected at elevations of 4600, 4800, 4950 and 5100m 
(Figure 7C). We collected samples at each location at 10:00 a.m. 
on following sunny days to avoid wide variations in weather and 
time. Also, we collected samples from within fences built in 2006 to 
be sure that samples had not been impacted by grazing animals. At 
each site, 100 to 200 g of healthy leaves from a 10 m^ region were 
randomly selected and immediately frozen by liquid nitrogen for 
later protein extraction and enzyme analysis. The samples were 
collected and the experiments were conducted in triplicate. 

Stomatal observation and measurement 

Environmental scanning electron microscopy (SEM) was 
conducted to evaluate the stomatal shape, size, and density using 
a QUANTA 200 3D (FEI, Tokyo, Japan) equipped with a large- 
field secondary electron detector and operated at 20 kV. Fresh 
specimens were frxed in formaldehyde: acetic acid: ethanol (FAA) 
at 4°C for 24 h. The samples were then dehydrated through a 
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Figure 7. Geographical location of the sampling site. A: The research site in the map of the eurasia continent. B: The geographical position of 
the research site on the Tibetan Plateau. C: Four sample sites selected at elevations of 4600, 4800, 4950 and 5100 m on the south-facing slope of the 
Nyainqentangiha Mountains. 
doi:1 0.1 371 /journal.pone.009841 0.g007 



graded ethanol series of 70, 85, 95, and 100% ethanol, each for 
45 min. Specimens were critical-point dried using liquid carbon 
dioxide and then coated with gold and palladium. With the SEM 
results, five leaves in the same general location were used to record 
stomatal density and the stomatal aperture length of 20 stomata 
was measured for each sample. 

Analysis of chlorophyll and carotenoid contents 

To measure the chlorophyll content, whole fresh leaves (1 g) 
were soaked in 10 ml acetoneiethanol (1:1) for 24 h under 
darkness. After centrifugation at 12 000 xg for 10 min, the 
supernatant was used for chlorophyll content measurement at 
absorbance of 645 and 663 nm using a spectrophotometer (Bio- 
Rad) [44]. 

Extraction and determination of malondialdehyde, 
proline, and abscisic acid content 

The malondialdehyde content was determined as described by 
Duan et al. [45]. Approximately 0.5 g of fresh leaves were 
homogenized in 10 mL of 10% trichloroacetic acid (TCA) and 
centrifuged at 12 000 for 10 min. Then, 2 mL of 0.6% 
thiobarbituric acid in 10% TCA were added to an aliquot of 2 mL 
of the supernatant. The mixture was heated in a boiling water bath 
for 30 min and then quickly cooled in an ice bath. After 
centrifugation at 1 0,000 x^ for 10 min, the absorbance of the 
supernatant at 450, 532, and 600 nm was determined. The MDA 
concentration can be estimated through the formula C (nmol 
ml~') = 6.45(A532-AgQQ) - O.56A450. The MDA concentration was 
expressed as |a.mol g ' fresh weight. 

The proline content was measured as previously reported [46] . 
Approximately 0.5 g of fresh leaves were homogenized in 8 mL of 
3% aqueous sulphosalicyhc acid, and the homogenate was 
centrifuged at 2 000 x^ for 1 0 min. Then, 2 mL of the extract 
were reacted with 2 mL of acidic-ninhydrine and 2 mL of glacial 



acetic acid for 1 h in boiling water. The reaction mixture was 
extracted with 4 mL toluene. The chromophore containing 
toluene was separated and the absorbency read at 520 nm with 
a spectrophotometer (Bio-Rad). 

The abscisic acid quantification was performed using a 
previously described method [47]. Lyophihzed samples from 
different regions were ground to a fine powder in liquid nitrogen 
using a mortar and pestle. Triplicate samples (50 mg dry weight 
each) were extracted with 5 mL of 80% acetone containing 
100 mg/mL 2, 6-ditert-butyl-methyl phenol and 500 mg/L citric 
acid for 16 h at 4°C in the dark. The extracts were further 
homogenized using a Polytron (Brinkmann Instruments, West- 
bury, NY, USA) at maximum speed for 1 min and centrifuged at 
3000 x.§ for 5 min. 

Antioxidant enzyme activity assays 

The activities of catalase, ascorbate peroxidase, glutathione 
reductase, and superoxide dismutase were measured spectropho- 
tometricaUy by monitoring the change in A240, A290 and A560, 
respectively [48,49]. A 5 g leaf sample was homogenized to a fine 
powder with a mortar and pestle under liquid nitrogen. Then, the 
powder and pestle were ground with 30 mL of extraction buffer 
containing 50 mM sodium phosphate buffer (pH 7.0), 0.2 mM 
EDTA and 2% polyvinylpolypyrrolidone (PVPP) for 10 min. The 
homogenates were filtered through two layers of cheesecloth and 
centrifuged at 4°C at 15 000 xg for 15 min. Supernatants were 
desalted on a Sephadex G-50 column (Pharmacia) and this final 
sample was used to measure the enzymatic antioxidant activities. 

Protein extraction and two-dimensional gel 
electrophoresis 

Protein extraction and two-dimensional electrophoresis were 
performed according to a previous method [50], with minor 
medications. Approximately 10 g of leaves from samples collected 
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at elevations of 4600, 4800, 4950, and 5100 m were ground in 
liquid nitrogen and total soluble proteins were extracted on ice in 
acetone containing 10% (w/v) TCA and 0.07% (w/v) dithiothre- 
itol (DTT). The liomogenates were held at — 20°C for 4 h and 
then centrifuged (8 000 x^, 30 min, 4°C). The pellets were washed 
with acetone containing 0.07% (w/v) DTT at — 20°C for 30 min 
and then centrifuged (8 000 x^, 20 min, 4°C) a total of three times. 
Finally the pellets were vacuum-dried and then dissolved in lysate 
(7 M urea, 2 M thiourea, 4% (w/v) CHAPS (3-[(3-cholamidopro- 
pyl)-dimethylammonio]-l -propane sulfonate), and 60 mM DTT) 
for 2 h at room temperature with intermittent shocking, and then 
the samples were centrifuged (12 000 x^, 20 min, 20°C). The 
supernatants were collected for 2-DE experiments with 900 |j.g of 
total proteins by a method used previously [51-53], which were 
executed in triplicate. 

In-gel digestion and MALDI-TOF/TOF analysis 

Protein spots that showed significant changes in expression with 
change in elevation were excised manually from colloidal CBB- 
stained 2-DE gels. Protein digestion with trypsin was first 

performed; then mass spectrometry analyses were conducted 
using a MALDI-TOF/TOF mass spectrometer 4800-plus Pro- 
teomics Analyzer (Applied Biosystems, Farmington, MA, USA) 
according to methods previously described [51-53]. 

Database search 

The primary and secondary MS data were transferred into 
Excel files and used as inputs to search against an NCBI non- 
redundant database; the search was restricted to viridiplantae 
(green plants) using the MASCOT search engine (www. 
matrixscience.com). The search parameters were established as 
follows: no restriction of protein molecular weight; one missed 
trypsin cleavage allowed; cysteine treated by iodoacetamide; and 
oxidation of methionine. The peptide tolerance was 100 ppm and 
the MS/MS tolerance was 0.25 kD. Protein identifications were 
validated manually, with at least four peptides matching. The 
keratin contamination was removed and the MOWSE score 
threshold was greater than 40 (P<0.05). Only significant hits were 

References 

1. Billings WD, Mooney HA (1968) The Ecology of Arctic and Alpine Plants. 
Biological Reviews ofthe Cambridge Philosophical Society 43: 481-529. 

2. Diaz HF, Grosjean M, Graumlich L (2003) Climate variability and change in 
high elevation regions: Past, present and future. Climatic Change 59: 1^. 

3. Gottfried M, Pauli H, Grabherr G (1998) Prediction of vegetation patterns at the 
Umits of plant life: A new view of the alpine-nival ecotone. Arctic and Alpine 
Research 30: 207-221. 

4. Xiong LM, Schumakcr KS, Zhu JK (2002) Cell signaling during cold, drought, 
and salt stress. Plant Cell 14: S165-S183. 

5. Ma L, Sun XD, Yang YQ, Yang SH, Kong XX, et al. (2013) Comparative 
proteome analyses reveal that the adaptation strategies of the alpine plant 
Potmtilla saundersima along the altitude gradient in the Northwestern Tibet 
Plateau. Unpublished. 

6. Gepstein S, Glick BR (2013) Strategies to ameliorate abiotic stress-induced plant 
senescence. Plant Molecular Biology 82: 623-633. 

7. Miehc G, Mlehe S, Kaiser K, LiuJQ, Zhao XQ,(2008) Status and dynamics of 
Kobresia pygmaea ecosystem on the Tibetan plateau. Ambio 37: 272-279. 

8. Miehe G, Miehe S, Bach K, NoUing J, Hanspach J, et al. (2011) Plant 
communities of central Tibetan pastures in the Alpine Steppe/ Kobresia pygmaea 
ecotone. Journal of Arid Environments 75: 711-723. 

9. Taylor SH, Franks PJ, Hulme SP, Spriggs E, Christin PA, et al. (2012) 
Photosynthetic pathway and ecological adaptation explain stomatal trait 
diversity amongst grasses. New Phytologist 193: 387-396. 

10. Beerling DJ, Chaloner WG (1993) Stomatal Density Responses of Egyptian 
Olea-Europaea L Leaves to Go2 Change since 1327 Be. Annals of Botany 71: 
431^35. 

1 1 . Fraser LH, Greenall A, Carlyle C, Turkington R, Friedman CR (2009) Adaptive 
phenotypie plasticity of Pseudoroegneria spicata: response of stomatal density, 
leaf area and biomass to changes in water supply and increased temperature. 
Annals of Botany 103: 769-775. 



accepted for the identification of the protein sample based on 
MASCOT probability analysis. 

Western blotting analysis 

SDS-PAGE was performed as described previously [54] using 
12% (w/v) polyacrylamide slab gels. For Western blot analysis, the 
protein samples were electroblotted onto polyvinylidene difluoride 
(PVDF) membranes using a Trans-Blot cell (Bio-Rad). After 
transfer, the membranes were probed with the appropriate 
primary antibodies and HRP-conjugated goat anti-rabbit second- 
ary antibody (Promega, Madison, WI, USA) and the signals were 
detected using an ECL kit (GE Company, Indianapolis, IN, USA). 
The primary antibodies (all obtained from Agrisera, Inc., Vannas, 
Sweden) were diluted as follows: polyclonal antibody against plant 
9-cis-epoxycarotenoid dioxygenase (NCED; 1:3000), dehydrin 
(1:3000), heat shock protein 18.2 (HSP18.2; 1:2000), mitogen- 
activated protein kinase 6 (MAPK6; 1:1000), and actin (1:2000). 

Statistical analysis 

Statistical analyses were performed using the statistical Software 
Package for Social Science (SPSS) version 12.0. One-way 
ANOVA for all variables was used for testing the treatment 
differences. Differences were considered significant at the P<0.05. 

Supporting Information 

FUe SI Supporting Figures. Figure SI, The second set 2-DE 
of four K. pygmaea samples from different elevations. Figure S2, 
The third set 2-DE of four K. pygmaea samples from difiFerent 
elevations. 
(ZIP) 

Author Contributions 

Conceived and designed the experiments: YPY XH. Perlbrmed the 
experiments: XL YQY LM XS SY XK. Analyzed the data: XL YQY. 
Contributed reagents/ materials/ analysis tools: YPY XH. Wrote the paper: 
XL YQY. 



12. Reddy KR, Robana RR, Hodges HF, Liu XJ, McKinionJM (1998) Interactions 
of C02 enrichment and temperature on cotton growth and leaf characteristics. 
Environmental and Experimental Botany 39: 117—129. 

13. Bettarini I, Vaccari FP, Miglictta F (1998) Elevated C02 concentrations and 
stomatal density: observations from 17 plant species growing in a G02 spring in 
central Italy. Global Change Biology 4: 17-22. 

14. Pcnuelas J, Matamala R (1990) Changes in N and S Leaf' Content, Stomatal 
Density and Specific Leaf- Area of 14 Plant-Spccics during the Last 3 Centuries 
of Co2 Increase. Journal of Experimental Botany 41: 1119—1124. 

15. GdXt J (1972) Availability of Carbon-Dioxide for Photosynthesis at High- 
Altitudes - Theoretical Considerations. Ecology 53: 494— &. 

16. Hultine KR, Marshall JD (2000) Altitude trends in conifer leaf morphology and 
stable carbon isotope composition. Oecologia 123: 32^0. 

17. Komer C, Bannister P, Mark AF (1986) Altitudinal Variation in Stomatal 
Conductance, Nitrogen-Content and Leaf Anatomy in Different Plant Life 
Forms in New-Zealand. Oecologia 69: 577-588. 

18. Qiang W'T, Wang XL, (.^hen 'L, Feng HY, An LZ, et al. (2003) Variations of 
stomatal density and carbon isotope values of Picea crassifolia at different 
altitudes in the Qilian Mountains. I'rces-Structurr and Function 17: 258-262. 

19. Franks PJ, Farquhar CD (2001) The effect of exogenous abscisic acid on 
stomatal development, stomatal mechanics, and leaf gas exchange in 
Tradescantia virginiana. Plant Physiology 125: 935-942. 

20. Israelsson M, Siegel RS, Young J, Hashimoto M, Iba K, et al. (2006) Guard ceU 
ABA and C02 signaling network updates and Ca2+ sensor priming hypothesis. 
Current Opinion in Plant Biology 9: 654-663. 

21. Bergmann DC, Lukowitz W, Somerville CR (2004) Stomatal development and 
pattern controlled by a MAPKK kinase. Science 304: 1494-1497. 

22. Wang HC, Ngwenyama N, liu YD, Walker JC, Zhang SQ_ (2007) Stomatal 
development and patterning are regulated by environmentally responsive 
mitogen-activated protein kinases in Arabidopsis. Plant CeU 19: 63-73. 



PLCS ONE I www.plosone.org 



12 



June 2014 I Volume 9 | Issue 6 | e98410 



Kobresia pygmaea and Elevation in Central Tibet 



23. Akashi K, Yoshimura K, Nanasato Y, Takahara K, Munckagc Y, ct al. (2008) 
Wild plant resources for studying molecular mechanisms of" drought/ strong light 
stress tolerance. Plant Biotechnology 25: 257-263. 

24. Joshi PN, Pradhan MK, Biswal B (2008) Interaction of Light Absorbed by 
Phytochrome and UV-B Radiation in Modulating the Composition and 
Function of Carotenoid in Photosynthetic Apparatus of Wheat Leaves: Role 
of UV-A Photoreceptor. Research Journal of Biotechnology: 319-324. 

25. Salgucro A, Leon R, Mariorti A, dc la Alorcna B. \^-ga JM. cl al. (2005) UV-A 
mediated induction of carotenoid accumulation in Dunalirlla bardawil with 
retention of cell viability. Applied MicrobiologN' and Biotechnology' 66: 506-51 1. 

26. Apel K, Hirt H (2004) Reactive oxygen species: Metabolism, oxidative stress, 
and signal transduction. Annual Review of Plant Biology 55: 373—399. 

27. Mitder R (2002) Oxidative stress, antioxidants and stress tolerance. Trends in 
Plant Science 7: 405-410. 

2f5. Lee DH, Lee CB (2000) Chilling stress-induced changes of antioxidant enzymes 
in the leaves of cucumber: in gel enzyme activity assays. Plant Science 159: 75- 
85. 

29. Liu YF, Qi MF, Li TL (2012) Photosynthesis, photoinhibition, and antioxidant 
system in tomato leaves stressed by low night temperature and their subsequent 
recovery. Plant Science 196: 8-17. 

30. Lou B, Xu DD, Xu HX, Zhan W, Mao GM, et al. {201 1) Effect of high water 
temperature on growth, survival and antioxidant enzyme activities in the 
Japanese flounder Paralichthys olivaceus. African Journal of Agricultural 
Research 6: 2875-2882. 

31. Yu C, Huang SJ, Hu XM, Deng W, Xiong C, et al. (2013) Changes in 
photosynthesis, chlorophvU lluorcsccnce, and antioxidant enzymes of mulberry 
(Morus ssp.) in response to salinity and high-temperature stress. Biologia 68: 
404-413. 

32. Yang F, Wang Y, Miao LF (2010) Comparative physiological and proteomic 
responses to drought stress in two poplar species originating from different 
altitudes. Physiologia Plantarum 139: 388-^00. 

33. Maiini L, Gaston KJ, Prosser F, Hulme PE (2009) Contrasting response of 
native and alien plant species richness to environmental energy and human 
impact along alpine elevation gradients. (Global Ecology and Biogeography 18: 
652-661. 

34. Higashi K, Shiota H, Kamada H (1998) Patterns of expression of the genes for 
glutamine synthetase isoforms during somatic and zygotic embryogenesis in 
carrot. Plant and Cell Physiology 39: 418-424. 

35. Kuznetsov W, Shevyakova NI (1997) Stress responses of tobacco cells to high 
temperature and salinity. Proline accumulation and phosphorylation of 
polypeptides. Physiologia Plantarum 100: 320—326. 

36. Mihailova G, Petkova S, Buchel C, (ieorgieva K (2011) Desiccation of the 
resurrection plant Haberlea rhodopensis at high temperature. Photosynthesis 
Research 108: 5-13. 

37. Ashraf M, Foolad MR (2007) Roles of glycine betaine and proline in improving 
plant abiotic stress resistance. Envirormiental and Experimental Botany 59: 206- 
216. 

38. Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, et al. (2006) Crosstalk 
between abiotic and biotic stress responses: a current view from the points of 



convergence in the stress signaling networks. Current Opinion in Plant Biology 

9: 436^42. 

39. Kosova K, Vitamvas P, Prasil IT (2011) Expression of dehydrins in wheat and 
barley under different temperatures. Plant Science 180: 46-52. 

40. Siaussat D, Laparie M, Maria A, Renault D (20 1 3) Heat shock protein responses 
to salinity, food deprivation, and temperature in the invasive ground beetle 

Merizodus soledadinus at the Kerguelen Islands. Polar Biology 36: 201—209. 

41. Kato S, Yamagishi K, 'I'atsuzawa F, Suzuki K, 'I'akano S, et al. (1993) 
Identification of (.Cytoplasmic and Nuclear Low-Molecular- Weight Heat-Shock 
Proteins in Tomato Fruit. Plant and Cell Physiology 34: 367-370. 

42. Ichimura K, Shinozaki K, Lena G, Sheen J, Henry Y, et al. (2002) Mitogen- 
activated protein kinase cascades in plants: a new nomenclature. Trends in Plant 
Science 7: 301-308. 

43. Li X, Yang YQ, Ma L, Yang SH, Kong XX, et al. (2013) Protein levels analysis 

Kobresia pygmaea (Cyjieraceae) adaptation to diurnal enviromnent on Qinghai- 
I'ibrt Plateau. Unpublished. 

44. Lichienlhalcr HK (1987) Chlorophylls and C^arotrnoids - Pigments of 
Photosynthetic Biomembranes. Methods in Enzymology 148: 350—382. 

45. Duan BL, Lu YW. \"m CY, Junttila O, Li CY (2005) Physiological responses to 
drought and shade in two contrasting Picea asperata populations. Physiologia 
Plantarum 124: 476-484. 

46. Bates LS, Waldren RP, Teare ID (1973) Rapid Determination of Free Proline 
for Water-Stress Studies. Plant and Soil 39: 205-207. 

47. Zhang AY, Zhang J, Zhang JH, Ye NH, Zhang H, et al. (201 1) Nitric Oxide 
Mediates Brassinosteroid-Induced ABA Biosynthesis Involved in Oxidative 
Stress Tolerance in Maize Leaves. Plant and Cell Physiology 52: 181 192. 

48. Jiang MY, Zhang JH (2001) Eflbct of abscisic acid on active oxygen species, 
antioxidative defence system and oxidative damage in leaves of maize seedlings. 
Plant and CeU Physiology 42: 1265-1273. 

49. Nakano Y, Asada K (1981) Hydrogen-Peroxide Is Scavenged by Ascorbate- 
Specific Peroxidase in Spinach-Chloroplasts. Plant and Cell Physiology 22: 867- 
880. 

50. Damerval C, Devienne D, Zivy M, Ihiellement H (1986) Technical 
Improvements in Two-Dimensional Electrophoresis Increase the Level of 
Genetic-Variation Detected in Wheat-Seedling Proteins. Electrophoresis 7: 52— 
54. 

51. Bai X, Yang L, Yang Y, Ahmad P, Yang Y, et al. (2011) Deciphering die 
protective role of nitric oxide agsunst salt stress at the physiological and 
proteomic levels in maize. J Proteome Res 10: 4349-4364. 

52. Yang LM, Tian DO, Todd CD, Luo YM, Hu XY (2013) Comparative 
Proteome Analyses Reveal that Nitric Oxide Is an Important Signal Molecule in 
the Response of Rice to Aluminum Toxicity. Journal of Proteome Research 12: 
1316-1330. 

53. Yang Y, ChenJ, Liu Q, Ben C, Todd CD, et al. (2012) Comparative proteomic 
analysis of the thermotolcrant plant Portulaca oleracea acclimation to combined 
high temperature and humidity stress. J Proteome Res 11: 3605-3623. 

54. Laemmli UK, Beguin F, Gujerkel G (1970) A Factor Preventing Major Head 
Protein of Bacteriophage T4 from Random Ag^egation. Journal of Molecular 
Biology 47: 69-74. 



PLCS ONE I www.plosone.org 



13 



June 2014 I Volume 9 | Issue 6 | e98410 



